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Generation of Coefficients for prediction filter in encoder. 



The invention relates to a transmitter for transmitting a digital information 
signal via a transmission medium, including: 

- input means for receiving the digital information signal, 

- adaptive prediction filter means adapted to derive a prediction signal from the digital 
5 information signal in dependence of an array of prediction filter coefficients, 

- first signal combination means for combining the digital information signal and said 
prediction signal so as to obtain a residual signal, 

- encoding means encoding for encoding said residual signal so as to obtain an encoded g^gn^t, 

- coefficient generator means for generating an array of filter coefficients A[i] in response to 
10 the digital information signal, i being an integer for which holds 0 < i £ p and whereby p is a 

variable, 

- output means for supplying the encoded signal to an output terminal for transmission via a 
transmission medium. 

The invention further relates to a receiver for receiving a transmission signal 
15 and generating a digital information therefrom, and to a transmission method. 

A transmitting device and receiving device of the type defined in the opening 
paragraph is known from J. Audio Eng. Soc, Vol. 44, No. 9, pp. 706 - 719, September 1 996. 

20 The known transmitting device is intended for efficiently reducing the bit rate for the 
transmission of a digital information signal. Therefore, prior to encoding the digital 
information signal a predicted version of the digital information signal is subtracted from the 
digital information signal. The thus obtained residual signal is successively encoded in the 
encoder and transmitted via the transmission medium. The performance of the linear 

25 prediction filter is crucial for the coding gain of the encoder. The performance of the linear 
prediction filter is determined by the prediction filter coefficient. A common method of 
finding the prediction coefficients (a-parameters) is the auto correlation method. It appeared 
that a-parameters determined with the auto correlation method did not result in the optimum 
coding gain, 
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The invention aims at providing a transmitting device and receiving device 
having a more efficient method of transmitting and receiving the digital information. 
5 The transmitter in accordance with the invention is characterized in that the 

device further comprises 

- smoothing means for smoothing the array of coefficients A[i] so as to obtain the array of 
prediction coefficients for supplying to the adaptive prediction filter means. 

The receiver in accordance with the invention is characterized in that the 

10 receiver further comprises; 

* smoothing means for smoothing the array of coefficients A[i] so as to obtain the array of 
prediction coefficients for supplying to die adaptive prediction filter means* 

The invention is based on the following recognition As long as transmission 
bandwidth and storage capacity is limited, there will be need to increase the bit reduction of 

1 5 digital signals. It appeared that the acquired prediction filter coefficients (a-parameters) 
obtained by for example the auto-correlation method, does not result in an optimal coding 
gain During extensive searches on the a-parameter values, it appeared that better coding 
gains could be achieved by slightly modified a-parameters. Coding gain improvements of 
around 4% could be achieved. However, the complexity needed in the transmitting device for 

20 such an extensive search would be enormous, and thus not practical* 

After investigation of the difference between the a-parameters before and after 
the extensive search, it appeared that the a-parameters are roughly the same (see Figure 8). 
When zooming in at the first few a-parameters the a-parameters obtained by the auto- 
correlation method show a rippling behavior while the a-parameters resulting from the 

25 extensive search are a smoothed version of it (see Figure 9). 

Instead of executing an extensive search on the a-parameters so as to obtain the 
prediction filter coefficients, the a-parameters acquired by the auto correlation method are 
post-processed by a smoothing function. The improvement of the coding gain coding gain 
obtained with the smoothing function is comparable with the coding gain obtained with 

30 extensive search. Figure 10 shows the results for a S-minute DSD fragment, where the result 
was even an improvement of 4.5%. Figure 1 1 shows that the performance is comparable of 
both methods (the extensive search would gain another 0.4% for this fragment). The increase 
of the complexity of the arrangement to implement the smoothing method is negligible 
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compared with the increase of the complexity of the arrangement when implementing the 
- extensive search method. 

These and other objects of the invention will become apparent from and 
5 elucidated further with reference to the embodiments described in the following figure 

description in which 

figure 1 shows an embodiment of an SACD coder, 
figure 2 shows the coding gain of successive frames, 
figure 3:shows the coding gain of successive frames and the corresponding 
10 value of max \aj\, 

figure 4 shows the coding gain as a function of max |a/|, 

figure 4a shows a detailed view of the coding gain as a function of max |a/j 

shown in figure 4, 

figure 5 and 5a shows the coding gain and impulse response as a function of the 
1 5 prediction order, 

figure 6 and 6 a shows the coding gain and impulse response as a function of 
the prediction order after applying the first two steps of the method in accordance to the 
invention, 

figure 7 shows the coding gain of the same successive frames in figure 2 after 
20 applying the method in accordance to the invention, 

figure 7a shows the value of max |o4 of successive frames after applying the 
method in accordance to the invention, 

figure 7b shows the coding gain as a function of the max \aj[ after applying the 
method in accordance to the invention, 
25 figure 8 shows an the arrays of prediction filter coefficient before and after 

extensive search, 

figure 9 shows a detail of the arrays of prediction filter coefficients before and 
after the extensive search, 

figure 10 shows the coding gain obtained by using the first array of coefficients 
30 versus die coding gain obtained by using the smoothed version of the first array of 
coefficients, 

figure 1 1 shows the coding gain obtained by using the array of coefficients 
obtained by extensive search versus the coding gain obtained by using the smoothed version of 
the first array of coefficients, 
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figure 12 shows an embodiment of a transmitter in accordance with the 

invention, .and 

figure 1 3 shows an embodiment of a receiver in accordance with the invention. 

5 Introduction 

For "Super Audio CD" (SACD), the audio signals are lossless coded using framing, linear 
prediction and entropy coding. The coding gainj is heavily determined by the quality of 
prediction. It is therefore important to have an algorithm that optimally determines the 
10 prediction coefficients. Figure 1 shows a block diagram of a simplified scheme mat is used in 
a SACD encoder. 

Disclosed is a low-complexity method for generating prediction coefficients 

that result in an improved coding gain. 

The problem is that frames, for which high coding gain can be expected, often 
15 suffer from significantly low gains. Such a low coding gam on frames brings the average 

coding gain lower, consequently having adverse implication on the storage capacity. Figure 2 
shows the coding gain of successive frames illustrating the problem. In an excerpt when the 
coding gain is relatively high (average 2.7), we see clearly that there are frames with a very 

low coding gain (close to 1 .0). 
20 Interactively, selecting a lower prediction order would have avoided this 

problem. However, an important issue is that coding gain is only available once the entire 
coding process has taken place and the complexity of interactively/iteratively re-encoding the 
frame using lower order is high. 

25 Disclosed Solution 

Here, we disclose a detection mechanism, which is capable of signaling LPC 
coefficients rendering low coding gain, where higher gain was possible, prior to encoding of 
the entire frame. The detection is based on the dynamic range of the direct form prediction- 
filter coefficients, often referred to as a-parameters. Consider Figure 3, where the top trace is 

30 same as figure 2 and the lower trace shows the value of the maximum of the absolute values of 
the a-parameters . It is clear from this Figure that the occurrence of low coding gain has a 
correlation with the occurrence of high values of max N- Based on mis detection, a 3-stage 
process is disclosed to eliminate this problem. The result is prediction coefficients offering 
higher performance of the encoding process. 
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In Figure 4, we illustrate a slightly different version of Figure 2. Here, the 
. maximum of the absolute values of a-p^rameters, max is shown on the horizontal axis and 
the vertical axis shows the corresponding coding gain. The a-parameters are determined by 
5 using the auto-correlation method in combination with the Schur algorithm. A line is drawn at 
a value of 1 0. Notice that for cases, where the maximum of the absolute value of a-parameters 
is higher than 80 9 there isn't a single occurrence of a frame with high coding gain. 

If we zoom in to one of the "problem frames", i.e, a frame with low coding gain 
and look at the prediction gain as a function of the prediction order, then we see Figure 5. 

10 Figure 5 and Sa highlights two issues. Figure 5 shows the coding gain as a function of the 
prediction order. Figure 5a shows the impulse response for various prediction orders. One 
may conclude from figure 5 that a higher coding gain was possible, for instance at prediction 
order 80, the gain was maximum > Hie second issue is highlighted in Figure 5a. The impulse 
response only changes in scale but not in shape. It is interesting to interpret this in the spectral 

15 domain, if one realizes that the prediction gain is proportional to the non-flatness in the 

spectrum [Flanagan paper]. Thus all die flattening has already been achieved in lower order, 
when the order is increased, only the scale is changing while the spectral coloring (or the 
flattening in the inverse filtered domain) is not changing. This points to some ill conditioning. 
Next, we discuss the steps that we have taken to avoid this problem by inspecting the impulse 

20 response of die prediction filter, i.e. without having to go through the entire coding process. 
We observed that that there is a relation between low coding gains and the maximum absolute 
value of die a-parameters returned by the auto-correlation method. For frames with a good 
coding gain the value of the maximum absolute a-parameters is always below 10, and problem 
frames show values of 100 or even 600. 

25 

Step 1: Detection of Problem Frames 

The first step in the method in accordance to the invention is the detection of 
the problem frames, i.e., frames for which the coding gain is low, while a higher gain was 
30 possible. Accordingly, Step 1 includes the computation of the a-parameters based on the auto- 
correlation method. Next, a value T t representing the absolute maximum value of a, is 
computed on die basis of die a-parameters, A = {a )f a*..., a S2B }, such that, T = max \A\. A 
frame is declared as a problem frame if T is greater than a predefined value, for example 10.0. 
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Step 2: Noise Addition 

Some of these problems are eliminated if one adds a very low-level random- 
noise signal to the input prior to LPC analysis. Equivalently, the value of ll[0] 9 the first auto- 
correlation coefficient, can be raised by a very small amount. If die threshold of step 1 is 
5 exceeded, the auto-correlation function coefficient R[0) is modified such that 

*t0]=*[0]* (1.0 + 4.10^), 
and the a-parameters are recalculated in accordance with the new value of H[0]. 

Although such a minor modification solves the problem in most cases, there 
were frames for which the problem persisted, or even got worse (in terms of coding gain). 

10 

Step 3: Prediction-Order Reduction 

The Tnfl y i TrtllTn absolute value of the now acquired a-parameters is again a good 
detection mechanism to see if problems are to be expected. For these frames where the 
addition of a little noise doesn't solve the problem, another measure must be taken and that is 
IS to take a reduced order for LPC prediction. It is difficult, however, to select the optimal 
prediction order. A threshold that seems to achieve good results is: 

New Order = min ( (current order * 8 / 9) , 80 ) 

20 where the current order is the order of prediction with which the analysis has been done. 

In practice it appeared that with the combination of these three steps (detection 
via max(abs(a-parameters)), noise addition, and order reduction), the coding gain improved 
considerably for problem frames. More importantly, die good frames were not effected 
adversely with this mechanism. 

25 

Results 

Figure 6, shows the improvements that have been achieved after applying the first 2 steps of 
the procedure discussed earlier. This is in contrast with Figure 5 where the coding gain 
dropped sharply beyond SO^-order prediction. Now the a-parameters also exhibit a reasonable 
30 behavior in terms of the amplitude dynamics. As indicated earlier, even after applying the 2 
steps, some problem frames persist and the application of step 3 eliminates them as well, 

Finally, Figure 7, 7a and 7b, shows clearly that the bad frames have been completely 
eliminated. 
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We have discussed a three-stage approach to eliminate the bad frames entirely. 
.This has been achieved by conditioning the signal in such a way that they deliver better 
prediction coefficients. 

5 While, only Schur recursion is discussed in this paper, experiments show that 

identical problems exist when using the Gholesky decomposition. Also here, the disclosed 
solution helps to completely eliminate the ill-conditioned a-parameters. 

"Smoothing" of the a-parameters could deliver (about) the same results as the 
1 0 extensive search, but for a negligible complexity. 

A smoothing algorithm uses the a-parameters to calculate a new set of a- 
parameters. Possible implementations are could be done by for example FIR or IIR-fUtering. 
After some experiments a variant of the HR-method was used selected: 

IS /* Apar[0] .. Aparftpo-1] contain the a-parameters of the filter with prediction order po */ 
for(i = l;i<po-l;i-H-) 
{ 

Apar[i] - Apar[i-1] + 2*Apar[i] + Apar[i+1] 

} 

20 

As can be seen from this part of code, the algorithm is recursive because the 
value that is adapted is reused in the next calculation. It appeared that it worked better than the 

non-recursive approach. 

On average some 3% coding gain improvement can be reached when applying 

25 this smoothing to the a-parameters. Figure 10 shows die results for a 5-minute DSD fragment, 
where the result was even an improvement of 4.5%. The complexity of this method is 
negligible. As can be seen from Figure 1 0 there were no frames where smoothing decreased 
the coding gain. Comparable results were obtained for some other DSD-fragments, with the 
exception that for a very few high gain frames>( >3) the gain showed a minor decrease. 

30 A comparison is made between the coding gains resulting from the extensive 

search the ones resulting from smoothing of the a-parameters that came out of die auto 
correlation method. Figure 1 1 shows that the performance is comparable of both methods (tho 
extensive search would gain another 0.4% for this fragment To decrease the complexity of the 
arrangement it is obvious that the method using smoothing of the a-parameters is preferred. 
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10 



The disclosed method is also lower in complexity when compared to the 
alternative that the entire coding process has to be performed before the coding gain is known. 

Whilst the invention is described with reference to preferred embodiments 
thereof, it is to be understood that these are not limitative examples. Thus various 
rnodifications may become apparent to those skilled in the art, without departing from the 
scope of the invention, as defined by the claims. 

The word •comprising' does not exclude the presence of other elements or steps 
than those listed in a claim. Any reference signs do not limit the scope of the claims. The 
invention can be implemented by means of both hardware and software. Several "means may 
be represented by the same item of hardware. Further the invention lies in each and every 
novel feature or combination of features. 
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CLAIMS: 



1. 



A transmitting device for transmitting a digital information signal via a 



transmission medium, including: 

- input means for receiving the digital information signal, 

- adaptive prediction filter means adapted to derive a prediction signal from the digital 
5 information signal in dependence of an array of prediction filter coefficients. 



prediction signal so as to obtain a residual signal, 

- encoding means encoding for encoding said residual signal so as to obtain an encoded signal, 

- coefficient generator means for generating an array of filter coefficients A[i] in response to 
1 0 the digital information signal, i being an integer for which holds 0 < i < p and whereby p is a 

variable, 

• output means for supplying the encoded signal to an output terminal for transmission via a 
transmission medium, 

characterized in that the device further comprises 
IS - smoothing m^ng for smoothing the array of coefficients A[i] so as to obtain the array of 
prediction coefficients for supplying to the adaptive prediction filter means. 

2. Transmitter as claimed in claim 1, characterized in that the smoothing means 
comprises low pass filtering means for low pass filtering of the coefficients so as to obtain the 

20 coefficient signal. 

3. Transmitter as claimed in claim 2, characterized in that the low pass filtering 
means are in the form of a FIR filter. 

25 4. Transmitter as claimed in claim 2, characterized in that the low pass filtering 

means are in the form of an 0R filter. 

5. Transmitter as claimed in any of the preceding claim in the form of an 

arrangement for writing the encoded signal on a record carrier. 



- first signal combination means for combining the digital information signal and said 
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6. Method of transmitting a digital information signal via a transmission medium, 

comprising the following steps: 

- receiving the digital information signal, 

5 - deriving a prediction signal from the digital information signal in dependence of an array of 
prediction filter coefficients, 

- combining the digital information signal and said prediction signal so as to obtain a residual 
signal, 

- encoding said residual signal so as to obtain an encoded signal, 

10 - generating an array of coefficients A[i] in response to the digital information signal, i being 
an integer for which holds 0 <S i £ p, whereby p is a variable, 

- supplying the encoded signal to an output terminal for transmission via a transmission 
medium, 

characterized in that the method further comprises the step 
15 - smoothing the array of coefficients A[i] so as to obtain the array of prediction filter 
coefficients. 

7 Receiver for receiving a transmission signal and generating a digital 



20 



information therefrom, the receive comprising: 

- receiving means for receiving the transmission signal and retrieving an encoded signal 



therefrom, 

- decoding means for decoding the encoded signal so as to obtain a residual signal, 

- adaptive prediction filter means adapted to derive a prediction signal from the digital 
information signal in dependence of an array of prediction filter coefficients, 

25 - signal combination means for combining the residual signal and the prediction signal so as to 
obtain the digital information signal, 

- coefficient generator means for generating an array of filter coefficients A[i] in response to 
the digital information signal, i being an integer for which holds 0 < i Z p and whereby p is a 
variable, 

30 characterized in that the receiving device further comprises 

- smoothing means for smoothing the array of coefficients ACi] so as to obtain the array of 
prediction coefficients for supplying to the adaptive prediction filter means. 
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ABSTRACT: 



A transmitter is disclosed for transmitting a transmission signal via a 
transmission medium The transmitter derives a prediction signal from the digital information 
signal in dependence of an array of prediction filter coefficients. The array of prediction filter 
coefficients has been obtained by smoothing a first array of coefficients. Said first array of 
5 coefficients has been generated in response to the digital information signal. A residual signal 
has been obtained by combination of the digital information signal and the prediction signal. 
The residual signal is encoded so as to obtain an encoded signal. The encoded signal is 
transmitted via the transmission medium. 
(Fig. 1) 
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